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In models with extra dimensions, a black hole evaporates both in the bulk and on the visible 
brane, where standard model helds live. The exact emissivities of each particle species are needed 
to determine how the black hole decay proceeds. We compute and discuss the absorption cross- 
sections, the relative emissivities and the total power output of all known fields in the evaporation 
phase. Graviton emissivity is highly enhanced as the spacetime dimensionality increases. Therefore, 
a black hole loses a significant fraction of its mass in the bulk. This result has important consequences 
for the phenomenology of black holes in models with extra dimensions and black hole detection in 
particle colliders. 

PACS numbers: 


A black hole is a grey body with temperature propor¬ 
tional to its surface gravity P|. The black hole emis¬ 
sion spectrum depends crucially on the structure and 
dimensionality of the embedding spacetime. In mod¬ 
els with large or warped extra dimensions, stan¬ 
dard model fields, except the graviton, are constrained 
to propagate on a four-dimensional submanifold of the 
higher-dimensional spacetime. In this scenario, the rela¬ 
tive emissivities of the fields (greybody factors) are essen¬ 
tial to determine if the black hole evaporates mainly on 
the four-dimensional brane or in the higher-dimensional 
bulk. This is particularly important for models of low- 
energy scale gravity, where detection of subatomic black 
holes in particle colliders and ultrahigh-energy cosmic ray 
observatories is only possible if a consistent fraction of 
the initial black hole mass is channeled into brane fields 
0. Counting of brane vs. bulk degrees of freedom (d.o.f.) 
provides a naive argument in support of the dominance 
of brane over bulk emission 0 : Since all standard model 
fields carry a larger number of d.o.f. than the graviton, 
brane emission must dominate over bulk emission. How¬ 
ever, a large emissivity for bulk fields could invalidate 
this conclusion 0. If the probability of emitting spin-2 
quanta is much higher than the probability of emitting 
lower spin quanta, the black hole may evaporate mainly 
in the bulk. 

The relative emissivities per d.o.f. of a four¬ 
dimensional non-rotating black hole are 1, 0.3A 0.11 and 
0.01 for spin-0, -1/2, -1 and -2, respectively [7|. In that 
case, the graviton power loss is negligible compared to 
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the loss in other standard model channels. Since brane 
fields are constrained in four dimensions, the relative 
greybody factors (relative emissivities) for these fields 
are expected to approximately retain the above values 
in higher dimensions. The graviton emission is expected 
to be larger due to the increase in the number of its he- 
licity states. The helicity states of a massless particle in 
Z?-dimensions are given by the representation of the lit¬ 
tle group SO{D — 2), i.e. the group of spatial rotations 
preserving the particle direction of motion. For instance, 
a five-dimensional graviton has five helicity states, corre¬ 
sponding to the SO{3) group of orthogonal rotations to 
the direction of motion. 

Therefore, a conclusive statement on brane vs. bulk 
emission rates requires the knowledge of the greybody 
factors for all fields. The higher-dimensional emission 
rates for spin-0, -1/2 and -1 fields are known 0 . How¬ 
ever, the exact emission rate for spin-2 fields has not 
been computed yet. (See, however. Ref. 0 ). In this pa¬ 
per we fill this gap. We find that the graviton emissivity 
is highly enhanced as the spacetime dimensionality in¬ 
creases. Although this increase is not sufficient to lead 
to a domination of bulk emission over brane emission, a 
consistent fraction of the higher-dimensional black hole 
mass is lost in the bulk. This fact has important con¬ 
sequences for the phenomenology of black hole events in 
low-scale gravity and primordial black hole formation. 

The Z3-dimensional non-rotating black hole is de¬ 
scribed by the higher-dimensional Schwarzschild metric, 
also known as Tangherlini metric M 

ds'^ = -fdt'^4-f~^dr‘^-\-r‘^dnl)_2, ( 1 ) 

where is the line element on the unit sphere 5^“^ 

and the black hole radius has been set to = 1 without 
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loss of generality. The formalism to handle gravitational 
perturbations of the Tangherlini metric has been devel¬ 
oped by Kodama and Ishibashi [m following prev ious 
work in four dimensions by Regge and Wheeler and 
Zerilli 0. There are three kinds of gravitational per¬ 
turbations: scalar, vector and tensor. Scalar and vector 
perturbations have their counterparts in T) = 4. The 
tensor type appears only in higher dimensions. The evo¬ 
lution equation for the perturbations is 


by dr/dr^ = f. For vector and tensor perturbations, the 
potential V is 

Ujl + D-S) iD-2)iD-4) (l-p2)(j_2)2 

^2 4j,2 4^D-1 

( 3 ) 

where I > 2,3,... is the angular quantum number and 
p = 2 (0) for vector (tensor) perturbations, respectively. 
For scalar perturbations, the potential is 


^ + {u;^-V)^ = 0, (2) 

where r is a function of the tortoise coordinate r* defined 


V = 


r Q{r) 
16r2i7(r)2 ’ 


( 4 ) 


g(r) = {D- 2f{D - ifx^ + {D- 2){D - 1){4[2(D - 2f - 3(D -2) + 4]m + {D - 2){D - ‘i){D - 6){D - l)}x'^ 
- 12{D - 2)[{D - 6)m + {D - 2){D - 1){D - A)\mx + -f 4,0(0 - 2)m^ , (5) 

H{r) = m + ^{O — 2){0 — l)x , m = l{l + O — 3) — {O — 2), x = llr^~^. ( 6 ) 


Assuming a harmonic wave and ingoing waves near 
the horizon, we have the following boundary condition: 

4'(r) ^ - (r - ^ r^r+. (7) 

For large r* we have both out- and in-going waves: 

4'(r) ^ , r* ^ oo . (8) 

The absorption probability for the wave is |Ap = (|T|^ — 
i?P)/|rp. The number of absorbed particles is 


where $ is the flux and a is the absorption cross-section. 
A sum over all final states and an average over initial 
states is understood. Solving the equations for the grav¬ 
itational perturbations with an expansion in spherical 

I 


harmonics, we find 

a = C{0,uj) ^ ^Nis |Azs|^ + Niv \Aivf + Nit \AiTf‘ 

( 10 ) 

where the subscripts S, V and T refer to scalar, vector 
and tensor perturbations, respectively. The multiplici¬ 
ties Ni are defined in Ref. |l^ and the normalization is 
C{0,u}) = 2[2 tt/uj)^-'^[0{0 - 3)flD_2]“^ 

The wave equation, Eq. @, can be solved in the low- 
energy regime w ^ 1. (For details, see Ref. 0.) This 
method uses a matching procedure to find a solution valid 
throughout the whole spacetime for any value of p. Low 
frequencies give a substantial contribution to the total 
Hawking power emission, thus providing a good approx¬ 
imation of the exact result. The absorption probability 
is 


|A|^ = 1 - |7^P = 47r 


/WN D+21 — 2 

[2) 


F 


/, 2i+£(D^\^p A 21-p(D-2) Y 

2{D-3) j V 2{D-3) J 


( 11 ) 


The low-energy absorption probability for spin-0 fields 
|1(tI| and spin-2 tensor perturbations is recovered by set¬ 
ting p = 0 in the above equation. The result for vector 
perturbations is obtained by setting p = 2. Gravita¬ 


tional scalar perturbations cannot be dealt with analyt¬ 
ically. However, numerical simulations give an effective 
Pgravscaiar ^2 + 040~^. The low-energy absorp¬ 
tion cross section can be obtained from Eqs. Hill . For 
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instance, in four-dimensions, where the tensor contribu¬ 
tion disappears, the I = 2 mode gives the cross section 
ai ^2 = 47ra;'^/45. This result agrees with the well-known 
result of Ref. 0 • 

The wave equation can also be solved in the high- 
energy limit. In that case, the absorption cross section 
is expected to approximate the cross section for particle 
capture. This conjecture has been verified in a number 
of papers for the scalar field. (See Ref. 0 and refer¬ 
ences therein.) The proof for spin-1 and spin-2 fields 
is sketched below. (For further details, see Ref. M-) 

I 

dE dEg dEy dEx f dtu to 

dt dt dt dt J 27r 


where we have separated the individual contributions of 
each harmonics. The Hawking temperature is Tff = {D — 
3) / (47r). Note that the number of helicities is included 
in the multiplicity factors and the sum of scalar, vector 
and tensor contributions. 

The absorption probabilities A for all frequencies can 
be computed numerically. Equation m is integrated 
from a point near the horizon (typically r — 1 ~ 10“®), 
where the field behavior is given by Eq. O- The numer¬ 
ical result is compared to Eq. 0 at large r. A better 
accuracy is achieved by considering the next-to-leading 
order correction term (see Ref. M) 

^(r) ^ T (l + -) 6“*“’'* + i? (l - ^ oo . 

(13) 

This allows the determination of the coefficients T , R and 
the absorption probability A. The results for the total 
integrated power, Eq. (112), are summarized in Table U 
The values for the lower-spin fields living on the brane 
are taken from Ref. 0 . The graviton values are normal¬ 
ized to the four-dimensional case, where P = 1.52 x 10“®. 
This is in exact agreement with Page’s result. (See Ta¬ 
ble I of Ref. 0 .) The results for lower-spin fields are 
normalized to their four-dimensional values. It is worth 
discussing some features of the numerical results: (i) The 
relative contribution of the higher partial waves increases 
with D. For instance, in four dimensions the contribution 
of the I = 2 mode is two orders of magnitude larger than 
the contribution of the I = 3 mode. More energy is chan¬ 
neled in 1 = 3 mode than in ^ = 2 mode for D > 9. (The 
largest tensor contribution in ten dimensions comes from 
the I = 4 mode.) Contributions from high I are needed 
to obtain accurate results for large D. For instance, in 
ten dimensions the first 10 modes must be considered 
for a meaningful result. The values in Table Q] have a 
5% accuracy, (ii) The total power radiated in gravitons 


Since high frequencies can easily penetrate the gravita¬ 
tional potential barrier, the absorption probabilities of 
all fields approximate 1 as w ^ oo. The cross section in 
the high-energy limit must include the contribution from 
all I < uj. Therefore, the largest contribution to the cross 
section is given by high-Z modes. Since the high-Z limit of 
the wave equation is independent of the type of perturba¬ 
tion , it follows Ais = Aiv = Ait- The universality 
of the high-energy absorption cross section follows from 
Eq. n 01) and the properties of the multiplicities. 

The total energy flux for the gravitational radiation is 


-{NislAtff + Niv\A‘,ff + , ( 12 ) 


increases more rapidly than the power radiated in lower- 
spin fields as D increases. This is due to the increase 
in the multiplicity of the tensor perturbations, which is 
larger than the scalar multiplicity by a factor at high 
D. Therefore, the main contribution to the total power 
comes from the tensor (and vector) modes. Eor instance, 
in ten dimensions the tensor mode contributes roughly 
half of the total power output. 

Table II gives the fraction of radiated power per d.o.f. 
normalized to the scalar field, where the graviton value 
includes all the helicity states. In four dimensions, the 
power loss in gravitons is negligible compared to the 
power loss in lower-spin fields. The graviton channel is 
only about 5% of the scalar channel. This conclusion is 
reversed in higher dimensions. Eor instance, the gravi¬ 
ton loss is about 35 times higher than the scalar loss 
in Z? = 11. Although the graviton emission is highly 
enhanced, the large number of brane d.o.f. (more than 
100 for the standard model) assures that the brane chan¬ 
nel dominates on the bulk channel. However, power 
loss in the bulk is significant and cannot be neglected 
at high D] about 1/4 of the initial black hole mass is 
lost in the 11-dimensional bulk. The particle emission 
rates per d.o.f. are shown in Table HI. The relative emis¬ 
sion rates of different fields can be obtained by sum¬ 
ming on the brane d.o.f. Eor instance, the relative emis¬ 
sion rates of standard model charged leptons (12 d.o.f.) 
and the 11-dimensional bulk graviton are roughly 1:1. 
This ratio becomes ~ 40:1 in five dimensions. We find 
again that the bulk energy loss is significant in higher¬ 
dimensional spacetimes. This long overdue computa¬ 
tion is of paramount importance for the phenomenol¬ 
ogy of any higher-dimensional gravitational models. Our 
main result is that the power loss in the graviton chan¬ 
nel is highly enhanced in higher-dimensional spacetimes. 
This has important consequences for the detection of mi- 
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TABLE I: Total power P of Hawking radiation channeled into different fields. The first three rows correspond to fields 
propagating on the brane. The last row is the power radiated in bulk gravitons normalized to the four-dimensional case. 


D 

4 

5 

6 

7 

8 

9 

10 

11 

Scalars 

1 

8.94 

36 

99.8 

222 

429 

749 

1220 

Fermions 

1 

14.2 

59.5 

162 

352 

664 

1140 

1830 

Gauge Bosons 

1 

27.1 

144 

441 

1020 

2000 

3530 

5740 

Gravitons 

1 

103 

1036 

5121 

2 X lO"* 

7.1 X 10^ 

2.5 X 10® 

8 X 10® 


TABLE II: Fraction of radiated power per d.o.f. normalized 
to the scalar field. The graviton d.o.f. (number of helicity 
states) are included in the results. 


D 

4 

5 6 

7 

8 

9 

10 

11 

Scalars 

1 

1 1 

1 

1 

1 

1 

1 

Fermions 

0.55 

0.87 0.91 

0.89 

0.87 

0.85 

0.84 

0.82 

Gauge Bosons 

0.23 

0.69 0.91 

1.0 

1.04 

1.06 

1.06 

1.07 

Gravitons 

0.053 

0.61 1.5 

2.7 

4.8 

8.8 

17.7 

34.7 

TABLE HI: Fraction of emission rates per d.o.f. normalized 
to the scalar field. The graviton result includes all the helicity 
states and counts as one d.o.f. 

D 

4 

5 6 

7 

8 

9 

10 

11 

Scalars 

1 

1 1 

1 

1 

1 

1 

1 

Fermions 

0.37 

0.7 0.77 

0.78 

0.76 

0.74 

0.73 

0.71 

Gauge Bosons 

0.11 

0.45 0.69 

0.83 

0.91 

0.96 

0.99 

1.01 

Gravitons 

0.02 

0.2 0.6 

0.91 

1.9 

2.5 

5.1 

7.6 


croscopic black hole formation in particle colliders and 
ultrahigh-energy cosmic ray observatories, where a larger 
bulk emission implies larger missing energies and lower 
multiplicity in the visible channels. Despite the increase 
in graviton emissivity, a non-rotating black hole in the 
Schwarzschild phase will emit mostly on the brane due 
to the higher number of brane d.o.f. However, black hole 
energy loss in the bulk cannot be neglected in presence 
of extra dimensions. 

The effects of graviton emission may be even more 
dramatic when one includes rotation. Consider a rotat¬ 
ing four-dimensional black hole M- Graviton emission, 
which is supressed for small rotations, rapidly increases 
with the angular momentum J. (In four dimensions J 
ranges from 0 to M^). As J grows from 0 to 0.7 gravi¬ 
ton emissivity grows 3 orders of magnitude (see Table I 
in Ref. ) while the emissivities of fermion and gauge 
bosons grow less than one order of magnitude. A simi¬ 
lar behavior is expected in higher dimensions, where for 
D > 5 there is no upper bound on J, implying that gravi¬ 
ton emission dominates the evaporation process. Since 
known results for rotating black holes in D-dimensions 
do not include gravitons |2(|, this remains an open ques¬ 
tion. 
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